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ABSTRACT 

We present mid- and far-IR photometry of the high-redshift (z=4.69) dusty quasar 
BR1202-0725. The quasar was detected in the near-IR, at a flux level (0.7±0.2 mJy) 
consistent with an average Radio-Quiet Quasar at it's redshift. Only upper limits for 
the emission were obtained in the far-IR. These upper limits, when combined with data 
from ground-based telescopes, are the first direct evidence for a turn-over in the far-IR 
emission and hence confirm that a black-body dominates the SED at FIR wavelengths. 
This black-body is most probably cool dust, constrained to have a temperature below 
80K, for a /3 of 1.5. 
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1 INTRODUCTION 

It came as a surp rise w hen groups working with IRAS data, 
e.g. Soifer et al. ( [L984| ) & Lawrence et al. ( [l986| ), discovered 
galaxies with bolometric luminosities equivalent to those of 
quasars but emitting predominantly in the infra-red wave- 
length region. Even mo re su rprising was the discovery by 
Rowan-Robinson et al. (|l99l|) of IRAS FSC 10214-^4724, a 



galaxy at a redshift of z=2.29 with an apparent IR lumi- 
nosity of 4x10" Lq. While F10214-I-4724 has turned out to 



be gravitationally lensed (Eisenhardt et al. 1996), leading 



to a lower luminosity of 2x10^'^ Lq, it is still of extreme 
luminosity. 

Five high-redshift radio-quiet QSO's were selecte d, pri - 
marily from an APM catalogue, by McMahon et al. (199'1) 
because their similarity to F10214-f4724 in redshift and 
overall properties made them natural objects in which to 
look for IR-emission from dust. BR1202-0725 (see figure ^ 
for a DSS image of the quasar and surrounding field), at a 
redshift of 4.69, was the only one to be detected at 1.25 mm 
with IRAM. Assuming the detection was due to dust and 
that the dust had a sin gle temperature (80K was assumed), 
McMahon et al. (1994) derived a dust mass of the order of 
4x10* Mq , a considerable amount of dust to exist at such an 
early epoch. If the far-IR spectrum of BR1202-0725 is sim- 
ilar to that of F10214-I-4724 then it has an IR luminosity 
of the order of IO'^^Lq, making it one of the most luminous 
objects in the universe. 

Later IR observations of this object by Isaak et aL 



( p,994D at 450, 800 & 1100 iim, and by Benford et al. ( |199£| ) 
at 350 microns, supported a slightly lower dust temperature 
of 68K and 50±7 K respectively, although none of these ob- 



servations saw the turnover of the emission to confirm that 
it comes from a black-body. Indeed, its spectral index of 
2.7±0.4 between 1100 and 800/im, calculated by Isaak et al., 
does not formally rule out synchrotron-emission, although a 
spectral index above 2.5 is difficult to explain by anything 
other than cool dust. 

The detection of CO emission from approximately 
IxIO^^Mq of molecular hydrogen (comparable to th at in 
a present-day lumino us galaxy) by Ohta et al. (1996) and 
Omont et al. (1996b), and the detection by Omont et al. 



( 19961 ) of a close companion also containing large amounts 



of molecular gas, further reinforced the idea that the galaxy 
was undergoing a large burst of star formation, possibly in- 
teraction triggered. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observing Strategy 

The observations of BR1202-0725 had two goals: to try and 
detect the turnover in the far-IR emission; and to detect the 
near-IR emission from the continuum to better determine 
the underlying power-law. Early predictions from Isaak et 
al. (1994) were for flux densities in the 100 - 200/im range 



to reach of the order of 100 mJy, depending on dust temper- 
ature. We therefore proposed to make observations through 
the wideband C90 and C160 fibers (centred at 90 and 170 
/im) of the ISOPHOT instrument ( Lemke et al. 1996 ) on- 
board ESA's Infrared Space Observatory (ISO, Kessler et 
al. 1996). It was hoped these two far-IR observations would 
lie either side of the peak of the emission. 
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2.3.1 CAM Data Reduction 



Figure 1. The field around the quasar BR 1202-0725, taken from 
the DSS. The quasar is in the middle of the field and is indicated. 
The scale is consistent with figure ^ 



Table 1. ISO observations, taken on 14 July 1996. 

Observation TDT Wavelength Filter Exposure 

Mode number sees 



CAMOl 
PHT32 
PHT32 



24002102 8-15 

24002103 90 
24002101 170 



LWIO 348 
C90 606 
C160 2696 



The ISOCAM instrument, described in Cesarsky et al. 

(1996), was used to obtain a single observation through the 
LWIO (8 - 15 /im ) filter with the aim of detecting the near- 
IR emission of BR1202-0725. A normal RQQ quasar with 
the redshift of BR1202-0725 would have a flux of approx- 
imately ImJy in this band, and the ISOCAM observation 
was therefore planned to reach that sensitivity by employ- 
ing the microscanning mode. 



2.2 Observations 

Observations, listed in Table 0, were performed on July 14 
1996. The ISOCAM observation was taken in microscanning 
mode, with the field-of-view of the instrument scanned over 
the object of interest in a 3x3 raster pattern. The detector 
was read out every 2.1 seconds, leading to redundancy in 
the dataset. The redundancy is used to distinguish glitches 
and other artifacts from valid sources. The ISOPHOT ob- 
servations were taken as PHT32 ra sters of size 3x4. The 
ISOCAM (^iebenmorgen et al. 1999| ) and ISOPHOT ( |Lan 



reijs et al. 2000D Handbooks give complete descriptions of 
these observing modes. 



2.3 Data reduction 



The data was sliced and dark corrected in the standard way 
(Delaney 2000), making use of a time-dependent dark cor- 
rection. A thorough first-order deglitching was performed 
using an iterative sigma-clipping method (Metcalfe et al. in 
prep). 

No correction was made at this stage for the effects of re- 
sponsive transients as the existing correction algorithms can 
adversely affect the S/N for very faint sources. The transient 
correction was instead made at the end of the data reduction 
by applying an appropriate scaling factor. This factor was 
established by transient correcting reference sources identi- 
fied for this purpose in representative data sets. 

Dedicated faint-source processing follow ed essentially 
the method described in Altieri et al. ( 1998 ) and Metcalfe 
et al. (in prep). 

Long-term baseline drifts were removed by smoothing 
(via median filtering) the time-history of each individual de- 
tector pixel, as sampled by the hundreds of readouts made 
during the observation, and then subtracting the smoothed 
baseline from the nominal history. The images recorded at 
the several raster positions were organised into a cube so 
that all samples of a given sky-position line-up. For each 
sky position, this sky-position-vector is then sigma-clipped 
to give a second-order deglitching which has been found 
to be extremely effective in removing residual glitches. It 
should be noted that first-order deglitching acts on the time- 
history of a detector pixel, as it moves over the sky dur- 
ing a raster. Second-order deglitching acts on the history of 
samples of each sky position, and so takes implicit advan- 
tage of the sampling redundancy of the raster measurement 
technique. The two deglitching steps are therefore highly 
complementary. Following this deglitching stage the raster 
mosaic was constructed straightforwardly from the raster- 
position-images and standard aperture photometry could be 
performed on the resulting map. 

In order to calibrate the aperture photometry the data 
were re-processed after inserting a set of fake sources into the 
raw data cube. These consisted of theoretical model PSFs 
matched to the optical configuration in use. By perform- 
ing aperture photometry on these fake sources in a manner 
identical to that used for the actual source any flux-altering 
effects of the complex data-reduction algorithm could be 
calibrated. 

The fake-source insertion operation was repeated for 
a range of fake-source brightnesses. In this way the source 
strength in ADUgs (ADU per gain per second) in the raw 
data, corresponding to a given source strength in the re- 
duced map, was determined. 

At this point it was only necessary to scale the source 
strength in the raw data for the effects of responsive tran- 
sients (as described above) and to apply the standard cal- 
ibration scaling factor relating ADUgs to incident mjy for 
each CAM filter (Siebenmorgen et al. 1999), in order to ar- 
rive at an estimate of the flux of the target source. 



2.3.2 ISOPHOT Data Reduction 

The ISOPHOT P32 data was reduced using the PHOT In- 
teractive Analysis (PIA) package including a pre-release ver- 
sion of P32 processing routines developed by R. Tuffs (Tuffs 
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Figure 2. The ISOCAM image overplotted with contours from 
a DSS2 red image. The quasar is in the centre of the image, with 
the black hne above being the missing column 24. There is a slight 
shift between the two images of approximately 5". This offset is 
not significant as the uncertainty in the ISOCAM astrometry due 
to the wheel positioning jitter is greater than this. The position 
of the quasar is marked. 



et al. in prep.). Initial data processing (ERD to SCP) used 
the new P32 processing routines, with the later processing 
(SCP to AAP) using the standard processing routines. De- 
fault values for the reduction parameters were used in both 
the old and new routines. 

Maps were generated using the Trigrid interpolation 
and the first quartile normalization flat-fielding methods. 



3 RESULTS 

In the ISOCAM dataset, see Figure ||, BR1202-0725 was 
detected with a flux of 0.7 mjy. Deriving the scatter in pho- 
tometry of ten randomly placed apertures near the centre 
of the image produces a la value of 0.07 mJy, while the 
results of the simulations give a (random plus systematic) 
uncertainty of 0.2 mJy. We therefore take a value of 0.7±0.2 
mJy for the flux from this object in the 8-15 fj.m wavelength 
region. The background flux was determined to be 1.06±0.1 
mJy per sq.arcsec (45±4 MJy sr~^), consistent with the typ- 
ical background predictions for CAM LWIO given in the user 
manual. 

Other surrounding sources (visible on the DSS image) 
are also present in the ISOCAM image. As they are not in 
the centre of the image the observation was not optimised 
for them, leading to a higher uncertainty on their detections, 
and we therefore do not quote fluxes for them. 

No object was detected at the position of the quasar in 
either the 90//m or 170//m datasets. The ISOPHOT 170^im 
map is shown in figure ^. The effective noise in both maps 



Figure 3. The 170/im ISOPHOT map. The quasar is located in 
the centre of the image but has not been detected by ISOPHOT. 
Its position is marked. (A colour version of this image is available 
at the MNRAS website.) 



was derived by moving a box the size of one ClOO or C200 
array pixel over the two maps (avoiding the outer parts of 
the maps) and converting to Jy. This produced 1 sigma noise 
values of 16 and 15 mJy at 90 and 170/im - effectively the 
3 sigma upper limit for both maps is 50 mJy. 

Combined results of this and other papers are presented 
in Table ^. Figure ^ shows the IR data, along with a grey- 
body spectrum at a temperature of 68K. T his wa s derived 
using the method outlined in Benford et al. (1999), i.e. 



5"^ = B^Q,{1 - exp(-r)) 

with r = [u/vo)^ , where uq — 2.4 THz, the critical frequency 
at which the source becomes optically thin, and taking (3 = 
1.5. The current data does not support the temperature of 
50K, w ith the other parameters, presented in Benford et al. 
( 1999| ). It is, however, consistent with the higher tempera- 
ture of 68K as given in Isaak et al. (1994). The maximum 



temperature allowed by the 170/im upper limit is 80K, al- 
though then the fit to the longer wavelength data becomes 
poor. 

Figure 1^, shown as normalised against v, extends the 
previous figure to include the 6cm data. It can be compared 
with the average spectrum for a Radi o-Quiet Quasar given in 
either Figure 10 of Elvis et al. (1994) or Figure 7 of PoUetta 
et al. (2000). If we approximately normalise the fluxes of 
BR1202-0725 at 6cm and R to those of the average RQQ the 
15pm datapoint is also in agreement, whereas the 350 - 800 
pi,m fluxes from BR1202-0725 are higher than the average 
RQQ, indicating the presence of cool dust in this quasar. 
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Table 2. BR1202-0725 fluxes. 
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Figure 4. Data points from Table ^ plotted with a black-body 
of temperature 68K. The ISOCAM datapoint is plotted as a dia- 
mond, the ISOPHOT 3a- upper limits are plotted as squares and 
observations with other instruments are plotted as asterisks. The 
68K black-body is shown as a dash-dot line. 



Observed 
Wavelength 



Rest 

Wavelength 



Flux 



Precision 





/im, 


mJy 


mJy 




114 






11.5 


2.02 


0.7 


0.2 


90 


15.8 


< 50^ 




170 


30.6 


< 50^ 




350^ 


61.5 


106 


7 


450^ 


79.1 


92 


38 


800^ 


141 


50 


7 


llOO^ 


193 


21 


5 


1250^ 


220 


10.5 


1.5 


1250* 


220 


12.59 


2.28 


60000^ 


10500 


0.3 


0.2 




1. McMahon et al 

2. Benford et al 

3. Isaak et al. ( ^994|') 

4. Omont et al. (1996a|). 

5. These are 3cr upper limits. 

a level consistent with with originating in a normal quasar 
at the distance of BR1202-0725. 
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Figure 5. All data points from Table ^ plotted as relative lumi- 
nosity against frequency. The ISOCAM datapoint is plotted as a 
diamond, the ISOPHOT 3cr upper limits are plotted as squares 
and observations with other instruments are plotted as asterisks. 
The 68K black-body is shown as a dash -dot line, and the av- 
erage RQQ SED from Elvis et al. ( |l994| ) is shown as a dashed 
line. The IS/^m flux of BR1202-0725 is consistent with the aver- 
age RQQ, whereas it exhibits an excess of emission at the far-IR 
wavelengths. 



4 CONCLUSIONS 

We have detected the distant, dusty quasar BR1202-0725 
in the near-IR and obtained upper-limits for the flux at two 
far-IR wavelengths. The upper-limits for the flux levels at 90 
and 170/im, when combined with the previous ground-based 
measurements, are consistent with the far-IR and sub-mm 
being emitted from a grey-body, most probably arising from 
dust. The maximum temperature allowed by the 170/im up- 
per limit is 80 K, while the probable temperature is 68K. 
The near-IR flux level of 0.7±0.2 mJy at 11.5/im is at 
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